NASA  AVSCOM 

Technical  Memorandum  4054  Technical  Memorandum  88-B-013 


Predicted  and  Measured  Strain 
Responses  of  Isotropic 
Panels  to  Base  Excitation 


Karen  H.  Lyle 
Aerostructures  Directorate 
USAART A- AVSCOM 
Langley  Research  Center 
Hampton,  Virginia 

Jack  D.  Leatherwood  and  Edward  F.  Daniels 
Langley  Research  Center 
Hampton,  Virginia 


NASA 

National  Aeronautics 
and  Space  Administration 

Scientific  and  Technical 
information  Division 


1988 


supjii^ry 


The  objective  of  the  research  described  is  to  assess  the  accuracy  of  classical  linear  theory  for 
predicting  acceleration  and  strain  for  cantilevered  and  clamped^free-clainped-free  (O-F-C-F)  panels 
excited  through  the  base.  Aluminum,  stainless  steel,  and  titanium  panels  of  various  dimensions 
and  thicknesses  were  vibration  tested  by  using  a  broadband  random  signal  applied  through  a 
shaker  mounting  fixture.  The  strains  were  measured  at  nine  locations  on  each  cantilevered  panel 
and  at  five  locations  on  each  O-F-C-F  panel.  Accelerations  were  measured  on  the  base  and  at  one 
panel  location  for  both  the  cantilevered  and  the  C-F-O-F  panels.  Predictions  were  based  on  the 
Ritz  method  with  assumed  beam  functions.  The  measured  accelerations  of  the  base  and  modal 
damping  values  were  input  to  the  analysis.  Measured  and  predicted  modal  frequencies  agree  to 
within  -4  percent.  Comparisons  between  predicted  and  measured  strain  and  acceleration  spectra 
were  within  an  average  error  of  20  percent  for  both  the  cantilevered  and  the  C-F-O-F  panels,  \ 

-/y 

Introduction 


An  analytical  and  experimental  research  program  has  been  initiated  at  the  NASA  Langley 
Research  Center  to  investigate  the  strain  response  of  panels  under  severe  thermal  and/or  acoustic 
loads.  The  ultimate  goal  of  the  program  is  to  develop  methods  for  predicting  panel  stresses  and 
strains  that  may  lead  to  improved  understanding  of  sonic  fatigue  failure  in  advanced  aircraft  (high 
sound  levels  on  aluminum,  composite,  and  sandwich  panels)  and  in  advanced  aerospace  vehicles 
(high  temperature  and  high  sound  levels  on  sandwich  panels).  The  initial  research,  described 
herein,  focuses  on  flat,  rectangular,  isotropic  panels  at  ambient  temperatures. 

Many  acoustic-response  research  studies  have  been  conducted  on  conventional  aluminum  panels, 
and  the  resulting  large  base  of  experimental  data  has  enabled  good  empirical  understanding  and 
confident  sonic  fatigue  design  methods  (refs.  1  to  9).  More  recently  there  has  been  similar  research 
with  composite  and  sandwich  panels  to  develop  empirical  sonic  fatigue  methods  (refs.  10  to  14); 
however,  complications  in  predicting  the  structural  response  have  caused  the  design  methods  for 
composite  and  sandwich  panels  to  be  much  less  reliable  than  those  for  aluminum  panels.  These 
complications  include  poorly  predicted  natural  frequencies,  nonlinear  large  deflections,  and  frame 
and  panel  boundary  problems.  Because  sandwich  and  composite  panels  can  be  made  with  a  great 
variety  of  structural  properties,  a  reliable  analytical  model  is  highly  desirable  to  reduce  the  effort 
and  costs  associated  with  the  large  amount  of  test  data  needed  in  developing  empirical  design 
charts  similar  to  those  that  already  exist  for  aluminum  panels. 

Unfortunately,  poor  comparisons  between  measured  and  predicted  bending  and  membrane 
strains  have  occurred  in  the  composite-  and  sandwich-panel  research  programs  that  are  due  not 
only  to  the  structural  modeling  complexities  but  also  to  complications  in  noise  spatial  distribution 
and  frequency  spectrum.  These  complications  include  complex  incident  sound  waves  (i.e.,  neither 
plane  waves  nor  grazing  waves)  and  frequency  spectra  varying  greatly  over  the  frequency  range  of 
interest.  Although  the  predictions  of  natural  frequencies  and  deflections  compare  very  well  with 
measurements  and  the  predominance  of  the  fundamental  mode  can  be  established,  the  predictions 
of  bending  strain  can  differ  from  measurements  by  a  factor  of  2  or  more  (refs.  11  and  13). 

Because  of  these  past  difficulties  in  obtaining  agreement  between  measured  and  predicted 
strains,  a  program  assessing  the  accuracy  of  classical  linear  theory  for  predicting  strain  was  initiated 
in  a  controlled  environment.  The  dynamic  structural  excitation  was  restricted  to  low  levels  so  that 
the  use  of  small  deflection  theory  was  applicable  in  the  analysis. 

The  present  paper  presents  the  initial  results  of  this  work.  The  first  section  describes  the 
experimental  portion  of  the  research,  including  detailed  descriptions  of  the  panel  mounting 
fixture,  the  test  panels,  and  the  data  analysis.  The  next  section  describes  the  analytical  model 
which  applies  the  Ritz  method  with  base  excitation  to  cantilevered  and  clamped-free-clamped-free 
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(C-F-C-F)  panels.  These  sections  are  followed  by  the  comparison  of  experimental  and  analytical 
results.  The  final  section  summarizes  the  current  status  of  the  research. 
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base  acceleration,  in/sec2 
base  displacement,  in. 
eigenvector  for  A;th  eigenvalue 
panel  length,  in. 

fcth  complex  mode  shape  coefficient,  in. 
panel  width,  in. 

bending  stiffness,  Eh^ /[I2(l  —  /c2)],  in-lb 

Young's  modulus,  lb/in2 

displacement  function  for  panel 

frequency,  Hz 

panel  thickness,  in. 

panel  loading,  lb 

work  energy,  in-lb 

forced  excitation,  lb/in2 

kinetic  energy,  in-lb 

time,  sec 

displacement  of  base,  in. 
potential  energy,  in-lb 

panel  displacement  for  pseudouniform  pressure  excitation,  in. 

panel  displacement,  in. 

distance  along  panel  length,  in. 

distance  along  panel  width,  in. 

fcth  mode  shape 

coefficients  in  cantilevered/clamped-clamped  beam  function 

coefficients  in  free-free  beam  function 

strain  parallel  to  r-axis  (longitudinal),  in/in. 

strain  parallel  to  y-axis  (lateral),  in/in. 

percent  of  critical  damping 

fcth  eigenvalue 

Poisson's  ratio 

mass  density  lb  sec2/ir.4 

mth  cantilevered/clamped-clamped  mode  shape 
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Abbreviation: 

FFT  fast  Fourier  transform 

Notation: 

=  d2/dt2 

V  del  operator 

Partial  derivatives  with  respect  to  x  and  y  variables  are  shown  as  subscripts  following  a  comma; 
e.g.,  F.x  =  dF/dx. 

Description  of  Experiment 

The  responses  of  cantilevered  and  clamped-free-clamped-free  (C-F-C-F)  rectangular,  isotropic 
panels  subjected  to  base  excitation  applied  at  the  clamped  ends  were  experimentally  determined 
in  this  investigation.  A  photograph  showing  the  apparatus  used  to  study  panel  response  for 
cantilevered  panels  is  shown  in  figure  1.  Sketches  of  the  cantilevered  and  clamped-free-clamped-free 
(C-F-C-F)  test  apparatus  are  shown  in  figure  2.  Cantilevered  boundary  conditions  were  achieved 
by  placing  the  panels  between  two  clamping  blocks  as  shown  in  figures  1  and  2(a).  To  obtain  as 
near  a  clamped  condition  as  possible,  it  was  necessary  to  place  lead  inserts  between  the  upper  and 
lower  panel  surfaces  and  the  clamping  devices  (blocks).  These  lead  inserts,  approximately  1/8  in. 
thick,  provided  high  damping  at  the  clamped  end  and  tended  to  fill  voids  along  the  clamped 
boundary  surfaces.  The  upper  clamping  block  was  attached  to  the  lower  clamping  block  by  a 
series  of  bolts  extending  through  the  lead  inserts  and  test  specimen.  The  lower  clamping  block 
was  bolted  directly  to  the  moving  shaker  head.  Clamping  pressure  on  the  test  panels  was  generated 
by  applying  torque  to  the  bolts  connecting  the  upper  and  lower  clamping  blocks.  Considerable 
effort  was  devoted  to  obtaining  clamped  boundaries,  and  the  above  method  provided  the  most 
consistent  and  repeatable  results  for  the  cantilevered  panels. 

The  C-F-C-F  boundary  conditions  proved  to  be  more  difficult  to  obtain  experimentally  because 
attempts  to  clamp  opposite  ends  of  the  test  panels  introduced  mechanical  loads  and  prestresses  that 
modified  the  stiffness  characteristics  of  the  panels.  This  was  due  to  factors  such  as  misalignment 
of  the  opposing  clamping  devices,  initial  panel  curvatures,  and  local  panel  irregularities  at  the 
clamping  surfaces.  To  minimize  the  problem  of  mechanically  induced  preloads,  the  clamping 
method  illustrated  in  figure  2(b)  was  used.  This  method  utilized  an  aluminum  support  fixture 
containing  a  fixed  clamping  device  at  one  end  and  a  movable  clamping  device  at  the  other  end.  The 
movable  clamping  device  provided  vertical  and  rotational  adjustment  capability  for  alleviating  the 
preloads.  This  was  accomplished  by  observing  selected  panel  strain  gauge  outputs  while  performing 
clamping  adjustments  at  the  movable  end.  It  should  be  noted,  however,  that  it  was  not  possible 
to  obtain  conditions  for  which  the  panels  were  totally  free  of  preloads,  particularly  in  the  y-axis 
(lateral)  direction. 

Eight  cantilevered  and  three  C-F-C-F  thin,  metallic,  rectangular  panels  were  tested  in  this 
investigation.  The  panels,  made  of  aluminum,  stainless  steel,  and  titanium,  were  milled  from 
available  sheet  stock  and  no  attempt  was  made  to  impose  strict  dimensional  tolerances.  The 
material  properties,  measured  dimensions,  and  transducer  locations  for  each  panel  are  given  in 
table  I,  with  the  key  to  the  dimensions  given  in  figure  3.  For  each  cantilevered  panel,  total  surface 
strains  were  measured  at  nine  locations.  Only  data  obtained  at  locations  2,  3,  5,  and  6  are  presented 
in  'h's  paper.  Data  from  five  surface-strain  locations  were  obtained  for  the  C-F-C-F  panels,  but 
only  sample  data  at  location  1  are  presented.  Accelerations  were  measured  a*  one  location  on  each 
panel  as  v« c*!  as  at  the  base.  However,  only  acceleration  data  for  panel  2  are  presented.  The  data 
selected  for  presentation  are  typical  of  those  obtained  at  the  remaining  locations. 


A  diagram  of  the  instrumentation  setup  is  presented  in  figure  4.  A  selected  base  excitation 
was  applied  to  the  clamping  fixture  using  an  electromagnetic  shaker  having  a  rated  force  of  100  lb 
(peak  sine  wave),  a  dynamic  stroke  of  0.75  in.  (peak  to  peak),  and  a  rated  velocity  of  100  in/sec. 
The  useful  frequency  range  of  the  shaker  was  dc  (direct  current)  to  6500  Hz.  Random  base  (shaker 
head)  acceleration  having  approximately  uniform  spectral  density  over  the  frequency  range  from 
50  to  600  Hz  was  obtained  by  band  pass  filtering  and  equalization  of  the  output  of  a  random-noise 
generator.  Equalization  was  necessary  in  order  to  compensate  for  roll-off  of  shaker  performance 
with  increasing  frequency.  The  equalizer  output  was  applied  to  a  power  amplifier  whose  output 
voltage  signal  provided  the  drive  signal  to  the  shaker  armature.  A  sample  base  acceleration 
obtained  from  an  accelerometer  attached  to  the  upper  clamping  block  is  shown  in  figure  5  for 
cantilevered  panel  4.  (Acceleration  is  given  in  g  units  where  1  g  «  32.174  ft/sec^.) 

Analog  outputs  from  each  sensor  were  applied  to  a  digital  processing  system  for  spectrum 
and  time  series  analysis,  identification  of  modal  characteristics,  and  signal  storage.  Each  of  the 
measured  panel  response  signals  was  analyzed  to  determine  spectral  characteristics  of  panel  strain 
and  acceleration  response,  panel-acceleration  to  base-acceleration  transfer  functions,  and  panel 
strain  to  base-acceleration  transfer  functions.  A  standard  circle-fit  procedure  (ref.  15)  was  also 
applied  to  identify  the  modal  frequency  and  damping  characteristics  of  the  observed  response 
modes.  The  results  of  these  analyses  were  used  for  comparison  with  theoretical  narrow  band 
predictions  of  panel  strain  and  acceleration  response. 

Analysis 

The  prediction  model  for  both  the  cantilevered  and  the  C-F-C-F  panels  is  based  on  the  Ritz 
method.  In  each  case,  for  the  clamped  boundaries,  the  model  assumes  no  translational  or  rotational 
motion  of  the  panel  relative  to  the  clamping  surface;  and  for  the  free  boundaries,  the  moment  and 
shear  force  are  zero. 

For  an  isotropic,  rectangular  panel  the  equation  of  motion  is  given  as  (ref.  16) 

D  V4  w{x,  y ,  t)  +  ph  iv{x,  y,  t)  =  p{x,  y,  t )  (1) 

where  D  is  the  bending  stiffness,  p  is  the  panel  density,  h  is  the  panel  thickness,  p(x,  y,  t)  is  the 
external  panel  loading,  and  w(x,  y,  t)  is  the  panel  displacement.  In  the  case  of  a  moving  base  with 
no  external  loads,  then  p(x,  y,  t)  =  0.  For  harmonic  motion,  where 

w(x,y,t)  =  v>{x,y)  el2nft 


then  equation  (1)  becomes 


D  V4  w(x,  y)  -  4t2phf2  w(x ,  y)  =  0 


(2) 


where  /  is  the  frequency.  The  boundary  conditions  for  the  cantilevered  panel,  where  the  base 
displacement  is  Aw,  are  given  as 


w(0,y)  =  Aw 

w,x  (0,  y)  =  w,xx  (a,  y)  =  w,XXx  {a,  y)  =  0 

w,yy  (x,  0)  =  W.yyy  (x,  0)  =  W,yy  (x,  6)  =  W,yyy  (jE,  b)  =  0  , 


(3a) 


The  boundary  conditions  for  a  C-F-C-F  panel  are  the  same  as  those  for  a  cantilevered  panel  except 
along  the  x  =  a  edge,  where  the  conditions  are  described  by 


w{a,  y)  =  Aw  and  w,x{a,y)  =  0 


(3b) 
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Using  the  transformation  (ref.  17) 


w{x,y)  =  W(x,y)  +  U{x,y) 


(4) 


where  W{x,y)  and  U(x,y)  are  components  of  panel  displacement,  and  substituting  into  equa¬ 
tion  (2)  yields 


D  V4  W (z,  y)  -  Air^ph. /2  W{x,  y)  =  —D  V4  U (z,  y)  +  An^ph /2  U (z,  y) 


(5) 


The  displacement  W(x,y)  can  be  solved  for  classical  cantilevered  boundary  conditions  for  an 
excitation  given  by  the  right-hand  side  of  equation  (5).  The  panel  displacement  U(x,y)  must 
satisfy  the  following  conditions: 

W(0,y)  =  Aw-U(0,y) 

W,x  (0,  y)  =  0  -  U,x  (0,  y)  =  —  £7,^(0,  y) 

W^zz  (ai  2/)  =  0  —  U,xx  (a>  y)  —  —U,xx{ai  y ) 

W^IZZZ  (<b  y)  —  0  —  U,x XX  (o,  y)  =  —U, xzz(fl)  y) 

lU,yy  (z,0)  =  0  —  [7,yy(z,  0)  =  —  C/,yy  (z,  0) 

W,yyy  {x,  0)  =  0  —  V,yyy  (z,  0)  =  ~U,yyy{x,  0) 

W,yy  (z,  b)  -  0  -  U,yy  (z,  6)  =  -U, yy(z,  6) 
lU,yyy  (z,  6)  =  0  —  C/,yyy  (z,  6)  =  —U,yyy[x,  b )  J 


(6a) 


For  the  C-F-C-F  panels,  the  conditions  at  z  —  a  are 

W{a,y)  =  AW-  U(a,  y) 

W,x  (a,  y)  =  0  -  U,x  (a,  y)  =  -t/,x(a 

For  either  '■ase,  for  this  experiment  the  function  U(x,y)  is 

U  =  Constant  =  Aw 


,y) } 


(6b) 


(7) 


which  would  make  the  right-hand  sides  of  all  parts  of  equation  6  equal  to  zero.  Substitution  of  U 
into  equation  (5)  gives 


D  V4  W(x,y )  —  An^phf2  W(x,y)  =  47r2 phf2Aw 


(8) 


The  response  W(x,  y)  is  calculated  by  the  Ritz  method  which  is  based  on  minimizing  the 
following  conservation  of  energy  equation: 


T  -  (V  +Q)  =  Constant 


(9) 


where  T  is  the  kinetic  energy,  V  is  the  potential  energy,  and  Q  is  the  work  energy.  For  a  rectangular, 
isotropic  panel  the  potential  energy  in  bending  is  given  by  (ref.  18) 

v  =  T  Jo  L  (^F'xx^x,y^+F,yy^x,y^2~2^~^{F'xx^x'y^ ^(^^“[^(^y)]2})  dxdv 
°  °  (10) 
the  kinetic  energy  is  given  by 


T  =  2n2phf^[  [  F^{x,y)  dx 
JO  Jo 


dy 


(11) 
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and  the  work  energy  is  given  by 


where  F(x,  y)  is  the  displacement  function  for  the  panel  and  q(x ,  y )  is  the  panel  loading. 

The  forced  displacement  response  of  the  panel  is  given  by  the  following  summation  of  mode 
shapes: 

r 

W{x,y)  =  ^2  Bk  Zk(x,y)  (13) 

k=  1 

where  B ^  is  a  complex  coefficient,  and  for  each  mode  the  mode  shape  is  given  by 


P  9 

Zk(xiV)  =  EE  Amn  Pm{x)  V’n(jz) 

m— 1  n=l 


(14) 


The  displacement  functions  4>rn(x)  and  ipniv)  are  assumed  beam  modes  (ref.  18),  where  in  the 
z-direction 


<t>m(x)  =  cosh (Pmx/a)  -  cos (/3mx/a)  -  am[sinh(/?mz/a)  -  sin(/?mz/a)]  (15a) 


and  in  the  ^-direction 

V’l  =  1 

^2(y)  =  v/3[l-(2y/6)] 

rjjn{y)  =  cosh {nny/b)  +  cos(7 ny/b)  -  6n[smh{^ny/b)  +  sin(7nj//t)] 


(«  >  3)  J 


(15b) 


The  reader  should  note  that  the  different  boundary  conditions  for  the  cantilevered  and  the 
C-F-C-F  panels  in  the  z-direction  result  in  different  numerical  values  for  /?m  and  am  in  equa¬ 
tion  (15a).  The  values  for  am,  Pm,  f>n,  and  7 n  are  calculated  in  reference  19. 

The  eigenvalues  X ^  and  eigenvectors  A^-  are  calculated  from  the  free-vibration  response,  i.e., 

V 

Q  =  0  in  equation  (9).  Minimizing  the  conservation  of  energy  expression,  where  F(x,  y)  =  Z^x,  y), 
with  respect  to  A  -  results  in  the  standard  eigenvalue  problem  of  the  form 


([X]-A[I]){A}=0  (16) 

where  [I]  is  the  identity  matrix  and  the  matrix  [X]  and  the  vector  {A}  are  determined  from 

dT  dV 


dA$.  dA)- 

i]  i] 


=  0 


(17) 


The  solution  of  the  matrix  in  equation  (16)  for  the  eigenvalues  and  eigenvectors  will  lead  to  the 
natural  frequencies  and  mode  shapes.  Numerical  values  for  the  eigenvalues  and  eigenvectors  are 
given  in  reference  18.  The  natural  frequency  / ^  is  related  to  the  eigenvalue  by 
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Returning  to  the  forced  solution  and  minimizing  the  conservation  of  energy  expression,  where 
F(x,y)  =  W(x,  y),  with  respect  to  B k  yields 


\f  j  SBk 


dV  A  o  dT 
n  — ■  =  47n — — 
dB k  dBk 


The  loading  on  the  plate  is  given  by 


q{x,y)  =  4n2phf2A 


Thus,  the  work  energy  Q  is 


Q  =  4ir2phf2Ai 


nb 

W(x,y) 

j 


dx  dy 


Solving  equation  (19)  for  {B}  yields 


B’.  = 


2Awf2  {Afni&ml Pm) 

m—l 


ft  ~  f 2 


(undamped) 


where  Aw  is  frequency  dependent.  For  small  damping  values,  the  percent  of  critical  damping  <; 
can  be  included  by  the  addition  of  2 i$f  fk  to  the  denominator  of  equation  (22): 


2Awf2  2_,  {^rn\amlPm) 


B,= 


fi  +  ^ffk-f2 


(damped) 


where  i  =  y/—l.  Since  the  excitation  is  customarily  measured  as  base  acceleration  Aa.  where 
Aa  =  47r 2 f2Aw,  then 


Bk  = 


2 Aa  [Am^am/ 0m) 

m=l 


4n2(f?  +  2i<;ffk  -  f 2) 


(damped) 


Thus,  Bk  can  now  be  substituted  into  equation  (13). 

The  strain  is  calculated  from  the  displacement  in  the  z-direction,  by 


h  d2w(x,y) 


and  in  the  y-direction  by 


_  h  d2w(x,y) 
(y~~2  dy* 


The  acceleration  is  given  by 


y)  =  -47 r2/2  w(x,  y )  =  -47r2/2[lT (x,  y)  +  Aw]  (26) 

The  root-mean-square  (rms)  spectral  values  in  the  x-direction  are  defined  for  rms  excitation  to 

ex, rms  (/)  =  (fx4)1/2  (27) 

where  the  superscript  *  denotes  the  complex  conjugate.  Thus,  substitution  of  equation  (24)  into 
(27)  gives 


be 


fj.rms(/)  — 


h  \  d2 
2  1  <9x2 


1  2 


(Bk)rea.\  ^ k 

Ljfc=  1 


+ 


ZL 

<9x2 


(-^A:)imag 

L/c=l 


I2  4  V'2 


(28) 


where  the  subscripts  real  and  imag  refer  to  the  real  and  imaginary  parts,  respectively,  of  the 
quantity  in  parentheses.  Likewise,  it  can  be  shown  that 
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Finally,  the  overall  (denoted  by  the  subscript  oa)  rms  strains  and  accelerations  over  the 
frequency  range  from  25  to  500  Hz  were  calculated  by 


tx.oa  — 


fy.oa  — 


Woa  — 


r  500 

yi  fi,rms(/) 
•/= 25 

500 
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“rms  (/) 

■/= 25 
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(31) 


Results  and  Discussion 

For  each  panel,  the  measured  and  predicted  strain  spectra  were  compared  for  a  frequency 
resolution  of  1  Hz.  As  shown  in  figure  6,  the  measured  results  for  cantilevered  panel  4  verify  that 
the  panels  vibrated  in  the  linear-response  region  so  that  a  linear  analytical  model  was  justified. 
Each  panel  was  tested  at  three  excitation  levels  and  showed  similar  linearity  results.  The  highest 
excitation  level  was  chosen  for  the  comparisons,  since  it  provided  the  best  signal-to-noise  ratio  in 
the  experimental  results.  Measured  and  predicted  natural  frequencies  and  overall  rms  strains  are 
presented  in  figures  7  and  8.  Panel  geometry,  material  properties,  and  transducer  locations  are 
described  in  table  I.  Figure  3  shows  the  transducer  locations. 


The  comparisons  for  measured  and  predicted  natural  frequencies  are  presented  in  figure  7. 
Included  in  the  figure  is  a  table  of  the  measured  natural  frequencies  and  damping  values.  The 
damping  values  indicate  that  the  test  panels  were  lightly  damped.  The  unusually  high  damping 
value  for  the  second  mode  of  cantilevered  panel  3  has  not  been  explained;  however,  this  mode  is  still 
considered  to  be  lightly  damped.  The  agreement  between  the  measured  and  predicted  frequency 
values  is  within  5  percent,  with  the  cantilevered  and  C-F-C-F  panel  natural  frequencies  predicted 
equally  well.  The  small  discrepancies  may  be  attributed  to  variances  in  panel  material  properties 
and  the  inability  to  obtain  perfectly  clamped  edges. 

The  comparisons  between  experimental  and  analytical  overall  rms  strain  values  are  shown  in 
figure  8.  Overall  rms  strain  is  the  total  strain  at  a  point  due  to  the  integrated  effect  of  all  spectral 
strain  components,  as  calculated  in  equation  (8).  The  lines  on  the  figure  denote  the  ±50-percent 
and  ±  100-percent  accuracy  boundaries  where  accuracy  is  defined  as  the  percent  of  deviation  of 
predicted  values  from  the  measured  values.  It  is  observed  that  most  of  the  points  fall  within  the 
50-percent  boundaries  with  the  exception  of  the  data  for  cantilevered  panel  5  (denoted  by  ©)  and 
the  transverse  (y-axis)  strain  for  the  C-F-C-F  panel  (denoted  by  EB). 

The  source  of  error  in  the  predictions  of  overall  strain  level  for  panel  5  is  not  clear.  It  may, 
however,  be  related  to  a  possible  coupling  of  panel  response  with  the  base  motion.  This  coupling 
is  illustrated  in  figure  9  which  shows  a  sharp  peak  and  dip  in  base  acceleration  at  frequencies 
close,  but  not  identical,  to  the  fundamental  frequency  of  panel  5.  Use  of  this  base-acceleration 
spectrum  as  input  to  the  analytical  model  could  have  contributed  to  the  large  difference  between 
measured  and  predicted  strains.  The  scurce  of  this  coupling  could  not  be  determined  although 
it  did  occur  for  another  panel,  panel  7,  made  of  the  same  material  and,  like  panel  5,  was  one  of 
the  heavier  panels  tested.  The  discrepancies  between  measured  and  predicted  transverse  strains 
for  the  C-F-C-F  panels  probably  resulted  from  changes  in  panel  curvature  and/or  mechanical 
preloads  induced  by  the  adjustable  clamping  device.  Although  the  longitudinal  (z-axis)  strains 
could  be  approximately  nulled  by  adjustments  of  the  clamping  mechanism,  it  was  not  possible  to 
set  the  transverse  strains  to  zero  simultaneously.  In  fact,  generally  the  lateral  (y-axis)  static  strains 
significantly  increased  as  the  z  axis  strains  were  minimized  during  adjustment.  This  resulted  in 
preloads  that  likely  affected  panel  stiffness  properties  and  significantly  influenced  y-axis  strain 
response,  particularly  at  resonance. 

If  the  y-strain  measurements  for  the  C-F-C-F  panels  are  omitted  because  of  their  questionable 
accuracy,  then  the  data  of  figure  8  would  indicate  that  the  average  difference  between  measured 
and  predicted  overall  rms  strains  is  approximately  20  percent.  This  agreement  is  felt  to  be  very 
good,  particularly  in  light  of  the  uncertainties  in  experimentally  derived  damping  values  and  the 
inability  to  achieve  perfectly  clamped  boundary  conditions. 

Measured  and  predicted  strain  spectra  for  location  2  on  panel  4  are  shown  in  figure  10. 
Location  2  was  nearest  the  theoretical  highest  strain  location  on  the  cantilevered  panel,  i.e., 
perpendicular  to  the  middle  of  the  clamped  edge.  In  addition,  the  strains  in  the  y-direction  were 
measured  and  predicted,  as  shown  in  figure  11  for  location  5  on  panel  4.  The  y-direction  strain 
in  these  panels  was  the  lowest  since  the  panels  behave  much  like  cantilevered  beams  with  very 
little  bending  across  the  panel  width.  In  both  figures  10  and  11,  the  experimental  and  analytical 
strain  spectra  compare  well  in  both  magnitude  and  shape  over  the  entire  frequency  range.  The 
higher  mode  resonant  peaks  do  not  match  exactly,  but  this  deviation  is  to  be  expected  because  of 
the  added  stiffness  effect  of  the  Ritz  method.  All  modes  in  the  z-direction  are  excited  because  of 
unsymmetric  boundary  conditions,  i.e.,  clamped  and  free.  However,  the  base  motion  excites  the 
panel  like  a  uniform  pressure,  so  only  modes  of  odd  half-wavelengths  in  the  y-direction  appear. 
Thus,  the  second  panel  mode  is  not  excited.  The  differences  in  measured  and  predicted  panel 
responses  of  the  resonant  peaks  are  not  unexpected.  It  is  generally  difficult  to  predict  peaks 
accurately  that  are  proportional  to  the  damping  values  used  in  the  analysis.  The  slight  differences 
between  predicted  and  measured  modal  frequencies  for  the  fundamental  cantilevered  modes  are 
attributed  to  imperfections  in  the  test  panels  and  inexact  panel  clamping. 


Figure  12  shows  typical  predict'd  and  measured  strain  results  for  a  C-F-C-F  panel,  location  1 
on  panel  1.  In  this  case  the  model  also  predicts  the  magnitude  and  shape  of  the  measured  strain 
spectra  accurately.  Since  the  panel  has  nearly  symmetric  boundary  conditions  in  the  i-direction. 
only  the  first  mode  is  excited  in  the  frequency  range  investigated,  i.e..  from  25  to  500  Hz.  Any 
nonsymmetry  in  the  x-direction  is  due  to  panel  imperfections  and  inexact  clamping  at  the  panel 
edges.  The  additional  restraints  of  the  second  clamped  edge  cause  a  higher  fundamental  frequency 
for  the  C-F-C-F  panels  than  for  the  cantilevered  panels.  However,  like  the  cantilevered  panels,  the 
discrepancy  in  the  magnitude  of  the  response  at  the  fundamental  frequency  may  be  attributed  to 
inaccurate  modal  damping  measurements. 

The  panel  accelerations,  unlike  the  strains,  depend  directly  on  the  acceleration  of  the  base.  (See 
eq.  (19).)  Hence,  it  was  important  in  the  analysis  to  incorporate  base  effects  and  phase  information 
correctly.  A  typical  example  of  the  ability  of  the  analytical  model  to  predict  panel  acceleration 
response  is  shown  in  figure  13  for  panel  2.  It  is  seen  that  the  analysis  was  able  to  predict  both 
the  resonant  and  antiresonant  response  characteristics  of  the  panel.  These  results,  together  with 
the  strain-response  results,  illustrate  the  capability  of  the  model  to  predict  accurately  both  strain 
response  end  acceleration  response  at  any  specified  location  on  the  surface  of  an  isotropic  panel. 

Concluding  Remarks 

To  date,  much  work  has  been  done  in  evaluating  the  accuracy  of  a  linear  isotropic  response 
theory  for  plate  deflection  as  well  as  modal  frequencies  and  mode  shapes.  However,  little,  if  any, 
work  appears  to  have  assessed  the  accuracy,  or  even  applicability,  of  such  linear  theory  toward 
predicting  strains.  In  conclusion,  this  paper  shows  that  the  Ritz  method,  using  assumed  beam 
functions,  predicted  the  linear  strain  and  linear  acceleration  response  spectra  of  isotropic  panels 
with  an  accuracy  of  approximately  20  percent.  The  discrepancies  between  predicted  and  measured 
results  are  likely  due  to  uncertainties  in  the  experimentally  determined  damping  values  and  the 
inability  to  obtain  ideal  clamped  boundary  conditions.  In  addition,  the  analysis  was  effective 
in  predicting  the  detailed  spectrum  characteristics  of  panel  response  as  well  as  the  variations  in 
response  across  the  panel  surfaces.  Having  verified  the  linear-strain  prediction  theory,  the  research 
program  can  more  confidently  progress  to  predicting  more  complicated  effects,  such  as  nonlinear 
panel  response  and  elevated  temperatures  on  strain  response. 

NASA  Langley  Research  Center 
Hampton.  Virginia  23665-5225 
June  24.  1988 


References 

1.  Lassiter.  Leslie  W.:  and  Hess.  Robert  W.:  Calculated  and  Measured  Stresses  in  Simple  Panels  Subject  to 
Intense  Random  Acoustic  Loading  Including  the  Near  Noise  Field  of  a  Turbojet  Engine.  NACA  Rep.  1367.  1958. 
(Supersedes  NACA  TN  4076.) 

2.  Clarkson,  B.  L.:  and  Ford.  R.  D  :  The  Response  of  a  Typical  Ahcraft  Structure  to  Jet  Noise.  J.  Royal  Aeronaut. 
Soc..  vol.  66.  Jan.  1962,  pp.  31  40. 

3.  Douglas  Aircraft  Co..  Inc.:  Structural  Design  for  Acoustic  Fatigue.  ASD-TDR-63-820.  U.S.  Air  Force.  Oct.  1963. 

4.  Schneider.  C.  W.:  Development  of  an  Empirical  Method  for  Calculating  Dynamic  Stresses  in  Flat  Elastic  Panels. 
ER-9702.  Lockheed-Georgia  <  o..  Dec.  29,  1967. 

5.  Ballentine.  John  R.;  Rudder,  Fred  F  ..  Jr.:  Mathis.  Janies  T.:  and  Plumblee,  Harry  E..  Jr.:  Refinement  of  Sonic 
Fatigue  Structural  Design  Criteria.  AFFDL-TR-67-156.  U.S.  Air  Force.  Jan.  1968.  (Available  from  DTIC  at' 
AD  831  118.) 

6.  Clckson.  B.  L.:  Stresses  in  Skin  Panels  Subjected  to  Random  Acoustic  Loading.  Aeronaut.  J..  vol.  72.  no.  695. 
Nov.  1968.  pp.  1000  1010. 

7.  Thompson.  V  G.  R.:  and  Lambert.  R.  F’.:  Acoustic  Fatigue  Design  Data ,  Part  II.  AGARD-AG-162-Pt.  II. 
Nov.  1972. 


8.  Rudder.  F.  F..  Jr.;  and  Plumblee,  H.  E.,  Jr.:  Sonic  Fatigue  Design  Guide  for  Military  Aircraft.  AFFDL-TR-74-112. 
I’.S.  Air  Force,  May  1975.  (Available  from  DT1C  as  AD  B004  600L.) 

9.  Jacobson.  Marcus  J.;  Sonic  Fatigue  Design  Data  for  Bonded  Aluminum  Aircrujl  Structures.  AFFDL-TR-77-45, 
r.S.  Air  Force.  June  1977.  (Available  from  DTIC  as  AD  B020  118L.) 

10.  Jacobson,  Marcus  J.:  Advanced  Composite  Joints  Design  and  Acoustic  Fatigue  Characteristics.  AFFDI  -TR-71- 
126.  I’.S.  Air  Force.  July  1971.  (Available  from  DTIC  as  AD  890  037L.) 

11.  Wright.  G.  O  :  Dynamic  Behaviour  of  F-'bre  Reinforced  Plastic  Beams  and  Plates.  Ph.D.  Thesis,  Univ.  of 
Southampton.  May  1973. 

12.  Holehouse.  Ian:  Sonic  Fatigue  Design  Techniques  for  Advanced  Composite  Aircraft  Structures.  AFWAL-TR-80- 
3019.  U.S.  Air  Force.  Apr.  1980.  (Available  from  DTIC  as  AD  A090  553.) 

13.  Teh.  Chin  Eng:  Dynamic  Behaviour  and  Acoustic  Fatigue  of  Isotropic  and  Anisotropic  Panels  Under  Combined 
Acoustic  Excitation  and  Static  In-Plane  Compression.  Ph.D.  Thesis,  Univ.  of  Southampton,  Aug.  1981. 

14.  Holehouse.  Ian:  Sonic  Fatigue  Design  Techniques  for  Advanced  Composite  Airplane  Structures.  Ph.D.  Thesis. 
Univ.  of  Southampton.  1984. 

15.  Ewins.  D.  J.:  Modal  Testing:  Theory  and  Practice.  John  Wiley  &  Sons  Inc.,  c.1984. 

16.  Ashton.  J.  E.:  and  Whitney.  J.  M.:  Theory  of  Laminated  Plates.  Progress  in  Materials  Science  Series.  Volume  IV, 
Technomic  Publ.  Co..  Inc.,  c.1970. 

17.  Myint-U.  Tyn:  Partial  Differential  Equations  of  Mathematical  Physics,  Second  ed.  Elsevier  Science  Publ.  Co.,  Inc.. 
c.1980. 

18.  Young.  Dana:  Vibration  of  Rectangular  Plates  by  the  Ritz  Method.  J.  Appl.  Mech.,  vo).  17,  no.  4,  Dec.  1950, 
pp.  448  453. 

19.  Young.  Dana:  and  Felgar,  Robert  P..  Jr.:  Tables  of  Characteristic  Functions  Representing  Normal  Modes  of 
Vibration  of  a  Beam.  Publ.  No.  4913,  Eng.  Res.  Ser.  No.  44,  Bureau  of  Engineering  Research.  Univ.  of  Texas, 
July  1.  1949. 


O  cfi 

S  S 


O 


3 

-O 


a 


i 


a 

.2 

’55 

a 

a> 

S 


a> 

a 

2 

Oh 

3 

o 

■*-> 

eC 

2 

T3 

9 


I 


a 

I 

J3 


v 

o 

3 

T? 


,<U 

0) 


s 


o 

£} 

S 

>> 


CO 


bo 

cc 


-2 

<L> 


X3 

2 

> 

a 


s 

o 


■'T 

a 

i 

£ 

0.2617  x  10~a 
.4042 

.2617 

.8084 

.8084 

.4042 

.8084 

.2617 

-  3 

JD 

£ 

O 

X 

LOOiOppppiO 

OlOOOOiOOO 

-c  .a 

OTfNOlMTfOM 

inrrwwa’jmcj) 

oooooooo 

d 

92, 

in. 

p 

r*  d 
Zf>  .— < 

uO 

CO 

co  d 
V  .3 

TfHNNCpOOtCOO 

OOH^oOqoOO 

ciddddddci 

io  d 
<u  -3 

OOpCOCOppOOCO 

»-H  i-H  «-H  «“H 

•'T  d 
4)  .3 

COSNOP-^hO) 

OOiOCOiOiCWXiO 

o 

1?  .g 

OCNCNC^CMpOCS 

dciciddddd 

CM  fl 

<u  .3 

<M^hCV|OOOwCM 

OpppppOp 

S  .3 

O®iO0iOiflH0 

Ot'-t'-b'-t'-t'-Ot''- 

“a  e 

pooppppp 

dddciciddd 

<j  a 

QOOOOO-htT 

o  *  ’  ’  ’  -J  ’  ’ 

-c  .S 

pooppppp 

dddddddw 

<s  g 

NCONOOOWtW 

popppppo 

ddd^fTrod^ 

Material 

6  p  E  p  S 

B  H  3  0  S 

a  B  ,C  3  B 

1 ! !  i  ii !  ^  1 

<H-<oq<55H(5d< 

Panel 

HNiOTful^NOO 

• '  -■  j  *  •  ■-  ' 


\v. 


ft 


(b)  Clamped-free-clamped-free  apparatus. 
Figure  2.  Sketches  of  test  apparatus  used. 


Accelerometer 


Figure  4.  Diagram  of  instrumentation  setup. 


Figure  5.  Sample  base  acceleration  for  cantilevered  panel  4. 
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Figure  6.  Relationship  of  rms  microstrain  base  acceleration  for  cantilevered  panel  4. 


Figure  7.  Measured  and  predicted  natural  frequencies  for  cantilevered  and  C-F-C-F  panels. 
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Figure  10.  Measured  and  predicted  strain  spectra  for  iocation  2  on  cantilevered  panel  4. 


Figure  11.  Measured  and  predicted  strain  spectra  for  location  5  on  cantilevered  panel  4. 
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Figure  12.  Measured  and  predicted  strain  spectra  for  location  1  on  C-F-C-F  panel  2. 
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Figure  13.  Measured  and  predicted  acceleration  spectra  for  cantilevered  panel  2. 
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